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(54) Data converter with minimum phase fir filter and method for calculating filter coefficients. 

(57) A data converter employs at least one mini- 
mum phase FIR filter (132,160,230,240 or 252), 
and includes an analog-to-digital converter 
(220) for converting an incoming analog signal 
into a plurality of digital signal samples, fol- 
lowed by a minimum phase FIR filter to filter the 
digital signal samples. Alternatively, the data 
converter includes a digital-to-analog converter 
(170) preceded by a minimum phase FIR filter to 
filter a plurality of digital signal samples that are 
converted into an analog signal by the digital- 
to-analog converter. The data converter may 
include both an analog-to-digital converter 
(220) and a digital-to-analog converter' (170). In 
a preferred embodiment, the minimum phase 
FIR filter (132,160,230,240 or 252) is an optimum 
minimum phase FIR filter. A method for pre- 
cisely calculating the filter coefficients of an 
optimum minimum phase FIR filter is also dis- 
closed. 
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Technical Field 

This invention relates generally to echo canceling 
and particularly to minimum phase finite impulse re- 
sponse (FIR) filters for attenuating the echo signal in 
communications circuits and a method for calculating 
tap weights of a minimum phase FIR filter. 

Background of the Invention 

The phenomena of signal echo occurs when a 
signal passes through an interface from a first trans- 
mission medium to a second transmission medium. 
An impedance mismatch across the interface will 
cause a portion of the energy of the signal to be re- 
flected back toward the source of the signal as an 
echo signal. The remaining portion of the energy of 
the signal propagates along the second medium. 

In communication applications, it is desirable to 
eliminate the echo signal, for example, to improve the 
quality of audio signals or prevent the occurrence of 
errors in data transmission. Basic echo cancellation 
techniques are explained in: the Bell System Techni- 
cal Journal, Vol. 59, No. 2, February 1980, pp. 149- 
159. One echo elimination technique discussed in the 
above journal is to detect when a signal is being trans- 
mitted in one direction and to attenuate all signals 
passing in the other direction. Errors in detecting the 
presence or absence of a signal can cause signal clip- 
ping when this echo elimination technique is used. 
Another technique utilizes subtractive attenuation. 
To eliminate the echo, the echo is replicated, such as 
being synthesized by an adaptive filter, then subtract- 
ed from the signal. Subtractive attenuation is more 
signal transparent and reduces signal dipping ef- 
fects. 

Another form of subtractive attenuation is ach- 
ieved with an echo canceler having fixed coefficients 
designed to approximate the impulse response of the 
echo path. Since the impulse response of the echo 
path is not known, this technique results in imperfect 
echo canceling. More precise echo canceling is ach- 
ieved when the echo canceler is adaptive. In an adap- 
tive echo canceler, the tap weights of the echo can- 
celer filter are changed over time to replicate the im- 
pulse response of the unknown echo path. 

If a FIR filter structure is employed in an echo 
path, the filter would ideally require an infinite number 
of tap weights since the unknown echo path generally 
has an impulse response of infinite duration. An adap- 
tive infinite impulse response filter having only a few 
tap weights could be used to achieve the necessary 
echo cancellation. Adaptive infinite impulse response 
filters, however, may become unstable. As a compro- 
mise, an adaptive FIR filter, having a finite number of 
tap weights, has been employed in echo cancellation 
applications. Adaptive FIR filters do not perfectly 
cancel the echo signal. 



An adaptive echo canceler adapts to the impulse 
response of the echo path. When some of the ele- 
ments of the echo path are known, the known ele- 
ments can be individually optimized to shorten their 

5 impulse response length, that is, to reduce the num- 
ber of tap weights. The reduced number of tap 
weights of individual elements contribute to shorten 
the length of the echo canceler. A reduced number of 
tap weights in a digital filter results in a reduced num- 

10 ber of multiply and add functions to carry out the fil- 
tering. Some of the elements in the echo path are dig- 
ital filters. Digital filters introduce a group delay dis- 
tortion and a total delay. To optimize a filter, and 
thereby minimize its length, both the group delay dis- 

15 tort ion and the total delay must be minimized. 

A constant group delay distortion can be ach- 
ieved only by a linear phase filter. Small variations in 
group delay distortion can be compensated for in an 
adaptive equalizer. The group delay through a filter 

20 must be as small as possible in the pass band of the 
filter. Minimizing the group delay cannot be achieved 
in a linear phase filter. Since the impulse response of 
the echo path is generally infinite in duration, to per- 
fectly cancel an echo signal requires an echo cancel- 

25 er that has an infinite impulse response. A minimum 
phase FIR filter provides both the shortest filter 
length and the minimum group delay distortion in a 
single filter. A minimum phase FIR filter also has all 
of its zeros on or inside the unit circle in the complex 

30 z-plane. Equivalently, a minimum phase FIR filter has 
no zeros outside the unit circle in the complex z-pla- 
ne. 

To employ a minimum phase FIR filter in an echo 
path would provide enhanced echo cancellation in 

35 that the length of the echo canceler would be shorter 
than for prior art methods. A further improvement in 
echo cancellation can be achieved by employing in an 
echo path a filter designed in accordance with the 
minimum phase alternation theorem set forth in Chen 

40 and Parks, "Design of Optimal Minimum Phase FIR 
Filters," Signal Processing, Vol. 10, June 1986, pp. 
369-383. A minimum phase FIR filter designed to 
meet the minimum phase alternation theorem will 
hereafter be referred to as an optimum minimum 

45 phase FIR filter. An optimum minimum phase FIR fil- 
ter has the shortest filter length and minimum group 
delay for a given set of magnitude constraints on the 
ideal frequency response. Group delay distortion re- 
sults from nonlinear phase. Although an optimum 

so minimum phase FIR filter does not have perfectly lin- 
ear phase, small deviations within the passband are 
permitted and are compensated for in an adaptive 
equalizer. However, extensive use of minimum phase 
FIR filters has been prevented by computational dif- 

55 ficulties in precisely calculating the filter coefficients. 
Minimum phase FIR filters were a curiosity, interest- 
ing to theoreticians but of little practical application. 
There are several reports in the literature of at- 
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tempts to design minimum phase FIR filters. For ex- 
ample: Hermann and Schussler, "Design Of Non-re- 
cursive Digital Filters With Minimum Phase/' Elec- 
tronic Letters, 1970, Vol. 6, pp. 329-330; Kamp and 
Wellekens, "Optimal Design of Minimum Phase FIR 5 
Filters"; Chit and Mason, "Design of minimum phase 
FIR digital filters", IEE Proceedings, Vol. 135, Decem- 
ber 1988, pp. 258-264; and Parks and Burrus, Digital 
Filter Design , John Wiley and Sons, New York, New 
York, 1987. However, the skilled artisan can not pre- 10 
sently use these methods to design practical mini- 
mum phase FIR filters because of limitations that ren- 
der the methods unable to accommodate more com- 
plex filter designs of practical importance or not con- 
sistently reliable. 15 

Hermann and Schussler introduced a method for 
specifying a minimum phase FIR filter to convert a fil- 
ter specification into a linear phase filter and subse- 
quently the linear phase filter into a minimum phase 
FIR filter. The Hermann and Schussler method has at 20 
least two major shortcomings. Firstly, it is not suffi- 
ciently general; the method can be used to design 
two-band filters but cannot accommodate filters re- 
quiring more than one stop band. Secondly, infinite 
attenuation cannot be specified at specific frequen- 25 
cies. In other words, the exact location of a zero can- 
not be specified. Using the Hermann and Schussler 
method, the skilled artisan can specify only limited 
attenuation over a frequency band. 

Kamp and Wellekens reported a method for opti- 30 
mal design of minimum phase FIRf ilters. Aconstrain- 
ed approximation procedure was used to obtain the 
magnitude function and the transmission zeros in the 
stopband(s). The zeros inside the unit circle are cal- 
culated via a low-degree polynomial factorization. A 35 
straightforward exchange algorithm is presented 
which achieves the constrained approximation step. 
This method has the same shortcomings as the Her- 
mann and Schussler method. 

Chit and Mason described a method for designing 40 
minimum phase FIR digital filters. The filter coeffi- 
cients are determined through a least means squared 
(LMS) algorithm such that the realized filter coeffi- 
cients are produced at the final convergence. The 
cost functions in the adaptation of the LMS algorithm 45 
are minimized. These cost functions are the frequen- 
cy domain specifications of the desired filter. How- 
ever, convergence noise in the LMS calculations 
causes nonoptimal location of filter zeros. 

Parks and Burrus disclose a three step proce- so 
dure, similar to Hermann and Schussler, for designing 
minimum phase FIR digital filters. The Parks and Bur- 
rus method takes advantage of characteristics of the 
impulse response polynomial to reduce the order of a 
polynomial ultimately factored. For each zero cross- 55 
ing on the frequency axis, a double pair of complex 
zeros (four zeros) is located on the unit circle at that 
frequency. The four zeros form a fourth order polyno- 



mial which can be readily factored out of the impulse 
response polynomial, or a reduced order polynomial 
from which other double pairs of complex zeros have 
been factored. Factoring the remaining reduced order 
polynomial is relatively easier than factoring the im- 
pulse response polynomial. The order of the remain- 
ing polynomial to be factored will always be less than 
the order of the original impulse response. Narrow 
passband filters have nearly all of their zeros on the 
unit circle. These zeros could be readily factored out 
of the impulse response polynomial. Factoring out 
these zeros significantly reduces the order of the im- 
pulse response polynomial. 

Conversely, for narrow stopband filters, the im- 
pulse response is a very high order polynomial and 
few zeros are located on the unit circle. The Parks and 
Burrus method does not provide a significant reduc- 
tion in the order of the polynomial to be factored. Even 
when the zeros on the unit circle are factored out, the 
remaining polynomial is very close in order to the or- 
der of the impulse response polynomial, and the dif- 
ficulty of factoring the very high order polynomial is 
not overcome. 

It would be desirable to employ minimum phase 
FIR filters or optimum minimum phase FIR filters in 
an echo path to minimize the impulse response 
length and the group delay of the filter. A method for 
precisely calculating filter coefficients for a minimum 
phase FIR filter that overcomes the problems of the 
prior art would facilitate using such filters. Employing 
minimum phase FIR filters, or optimum minimum 
phase FIR filters, in an echo path contributes to min- 
imizing the length of an echo can eel er. The method 
would be sufficiently general to accommodate multi- 
ple stop bands and permit precise location of zeros. 

Summary of the Invention 

In accordance with an embodiment of the inven- 
tion, a data converter is disclosed employing at least 
one minimum phase FIR filter. The data converter in- 
cludes an analog-to-digital converter for converting 
an incoming analog signal into a plurality of digital sig- 
nal samples, followed by a minimum phase FIR filter 
to filter the digital signal samples. Alternatively, the 
data converter includes a digital-to-analog converter 
preceded by a minimum phase FIR filter to filter a 
plurality of digital signal samples that are converted 
into an analog signal by the digital-to-analog conver- 
ter. The data converter may include both analog-to- 
digital and digital-to-analog conversion. In a preferred 
embodiment, the minimum phase FIR filter is an op- 
timum minimum phase FIR filter. A method for pre- 
cisely calculating the filter coefficients of an optimum 
minimum phase FIR filter is also disclosed. 
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Brief Description of the Drawing 

FIG. 1 is a block diagram showing an environment 
in which an embodiment of the invention may be 
used; and 

FIG. 2 is a block diagram showing an environment 
in which an alternate embodiment of the inven- 
tion may be used. 

Detailed Description 

In Fig. 1, modem 110 is a full duplex data trans- 
ceiver system. A first digital bit stream, or first data 
stream, is sent through host interface 100 to transmit- 
ter 130 and is filtered by transmitter filter 132 to re- 
move out of band energy from the first data stream. 
In the transmitter, the first data stream is modulated 
by an appropriate method, relevant to the modem 
standard that is employed, such as frequency shift 
keying, quadrature amplitude modulation, or quadra- 
ture phase shift keying, to provide the transmitter out- 
put first data stream. Alternatively, the transmit sym- 
bols may be encoded by baseband techniques such 
as 2B1Q and 3B2T. The transmitter output first data 
stream is the input to digital low pass filter 160 and a 
first input to echo canceler 140. Filter 160 filters out 
low frequency noise from the first data stream and 
provides the filtered first data stream as input to dig- 
ital-to-analog converter 170. Digital-to-analog con- 
verter 1 70 converts the first data stream to an analog 
signal. The analog signal is filtered by analog recon- 
struction filter 180 and then conducted to hybrid 190 
and onto a communication network. 

A portion of the analog signal output from filter 
180 is reflected back toward host interface 100 and 
received through hybrid 190 along with an analog sig- 
nal from a remote modem. This combined received 
analog signal is filtered by analog antialiasing filter 
210 and then converted by analog-to-digital converter 
220 to a second digital bit stream or second data 
stream. The second data stream is filtered by digital 
low pass filter 230 and digital high pass filter 240 to 
remove the DC and power supply frequencies. The 
filtered second data stream output from filter 240 is 
summed with the output of echo canceler 140 in sum- 
ming junction 260. The output data stream from echo 
canceler 140, when summed in summing junction 260 
with the second data stream output from filter 240, 
cancels the reflected portion of the transmitted first 
data stream. The second data stream output from 
summing junction 260 provides the input to receiver 
250 and a second input to echo canceler 140. Receiv- 
er 250 demodulates the encoded information in the 
second data stream to recover the data stream trans- 
mitted by a remote modem. Receiver filter 252 re- 
moves out of band energy. 

In a preferred embodiment, one or more of the 
transmitter filter 132, digital lowpass filter 160, digital 



low-pass filter 230, digital highpass filter 240, and re- 
ceiver filter 252, in any combination, is or are a mini- 
mum phase FIR filter. In a further preferred embodi- 
ment, all of the filters are minimum phase FIR filters. 

5 In yet a further preferred embodiment, one or more of 
the filters in any combination are optimum minimum 
phase FIR filters. Generally receiver filter 252 and 
transmitter filter 132 are not simultaneously minimum 
phase FIR filters, or simultaneously optimum mini- 

w mum phase FIR filters. 

The invention will be described in greater detail 
with reference to the block diagram of FIG. 1 in which 
an embodiment of the invention may be used in a mo- 
dem connected to a switched network where the ana- 

15 log-to-digital conversion samples at a rate at or above 
the Nyquistrate. Most interface 100 is coupled to mo- 
dem 110. More specifically, host interface 100 is cou- 
pled to transmitter 130, which includes transmitter fil- 
ter 132. Transmitter 130 is coupled both to echo can- 

20 celer 140 and the digital-to-analog section 1 50 of mo- 
dem 110. The digital-to-analog section 150 is com- 
prised of series-coupled digital low-pass filter 160, 
digital-to-analog converter (DAC) 170, and analog re- 
construction filter 180. Collectively, transmitter 130 

25 and the digital-to-analog section 150 comprise trans- 
mit path 300. Filter 180 is coupled to hybrid 190 which 
is isolated from the switched network by isolation 
mechanism 192. A modem (not shown) with which 
modem 110 is communicating is coupled to a remote 

30 end of the network. Also coupled to hybrid 1 90 to re- 
ceive data transmitted by the remote end modem is 
analog-to-digital section 200. The analog-to-digital 
section is comprised of serially coupled analog antia- 
liasing filter 210, analog-to-digital converter 220, dig- 

35 ital low- pass filter 230 and a digital high- pass filter 
240. A coder-decoder employing the invention is not 
limited to modem applications, but may be used in 
such applications as telephone switching equipment, 
speech coding systems, and audio systems. Acoder- 

40 decoder may not include digital high pass filter 240. 
Additional applications, requiring an echo canceler, in 
which the invention may be used include but are not 
limited to public switched network modems, private 
line two-wire modems, Integrated Services Digital 

45 Network (ISDN) 2B+D subscriber loops, High-speed 
Digital Subscriber Line (HDSL) Asymmetric Digital 
Subscriber Line (ADSL), and Fiber Data Distributed 
Interfaces(FDDI). Additional applications not neces- 
sarily requiring echo cancellation that require digital 

so filtering to provide low group delay, include but are 
not limited to telephone line cards, speech and audio 
coding systems and more generally audio systems. 

The high-pass filter is coupled both to echo can- 
celer 140 and to a receiver 250. The output of digital 

55 high-pass filter 240 is summed with the output of echo 
canceler 140 by summing junction 260. The output of 
summing junction 260 is coupled as an input to both 
echo canceler 140 and receiver 250. Receiver 250 in- 
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eludes filter 252. Receiver 250 is coupled to host in- 
terface 100 as an input. Collectively, analog-to-digital 
section 200 and receiver 250 comprise receive path 
302. The digital-to-analog section 150 and the ana- 
log-to-digital section 200 form a coder-decoder 400, 5 
which is the front end of modem 110. 

Host interface 100 might consist of a data termi- 
nal, such as an RS232 data terminal, or a modem em- 
bedded in a computer mother board application such 
that the remote end modem interfaces directly with 10 
the host's microprocessor bus. The outgoing data 
from host interface 100 is transmitted by transmitter 
1 30 and is filtered by transmitter filter 1 32 to remove 
out-of-band energy from the modulated digital sam- 
ples. The function of filter 132 may be embedded in 15 
filter 160, or the function of filter 160 may be embed- 
ded in filter 132. When the function of one of the fil- 
ters is embedded in the other filter, the combined fil- 
ter is generally a minimum phase FIR filter; simulta- 
neously filter 252 is not a minimum phase FIR filter. 20 
The data from transmitter 130 are modulated into 
words, such as, by the well-known "quadrature-ampli- 
tude phase-shift keying" (QPSK) technique. The 
transmitter output is in digital form and is input to both 
echo canceler 140 and the digital-to-analog section 25 
150. 

The primary purpose of the digital low-pass filter 
160 is to increase the sampling rate such that the 
sampling rate of the DAC 170 is relatively high. The 
transmitter 1 30 operates at a relatively low rate, such 30 
as the Nyquist rate of the highest frequency of inter- 
est. The digital low-pass filter 160 is an interpolation 
filter that increases the sampling rate by a multiple, 
such as L, from the Nyquist sampling rate to L times 
the Nyquist sampling rate. Filter 160 is not required 35 
when the digital-to-analog section 150 is operating at 
the Nyquist rate. 

Analog reconstruction filter 180 smooths the 
step-wise continuous analog output of DAC 170 to 
provide a filtered analog signal to data access ar- 40 
rangement 190 for transmission over the switched 
network to a far-end modem. Data access arrange- 
ment 190 includes an isolation device 192 to provide 
high voltage isolation. Atransformer is shown but the 
isolation device need not be a transformer. Other iso- 45 
lation devices such as linear opto-couplers would suf- 
fice. 

Some of the energy transmitted by modem 110 to 
data access arrangement 1 90 and the remote modem 
is reflected back into receiver 250. The remote mo- so 
dem is also simultaneously sending data that is re- 
ceived by modem 110 receiver 250. Some of the en- 
ergy transmitted by the remote modem is reflected 
back toward the remote modem as a reflected signal. 
The remainder of the energy that enters modem 110 55 
while transmitter 130 is active should be canceled by 
echo canceler 140. Ideally, the output of summing 
junction 260 contains only the signal from the remote 



end and not energy from transmitter 130 which is re- 
flected through data access arrangement 1 90. That 
is, the output of summing junction 260 does not con- 
tain a reflected component The transmitter 130 out- 
put that is reflected back into modem 110 is due to a 
mismatch of impedances, due to both the two-to-four 
wire conversion and the unknown load impedance of 
the line, a mismatch of unknown magnitude that nec- 
essitates the echo canceler. 

The analog antialiasing filter 21 0 is a digital filter 
intended to remove noise near the sampling rate fre- 
quency of the analog-to-digital converter 220. When 
the analog-to-digital converter 220 is sampling the re- 
ceived signal at a relatively high rate, typically some 
multiple of the Nyquist rate, the purpose of the digital 
low- pass filter 230 is to lower the sampling rate down 
to the Nyquist rate without introducing noise. The nar- 
row-band antialiasing filter 230 may have several sta- 
ges to achieve the filtering. Filter 230 is not required 
when the analog-to-digital section 200 is operating at 
the Nyquist rate. 

Each time a stage of filtering occurs, out of band 
energy is removed. Upon increasing the sample rate, 
the out of band energy is removed subsequent to in- 
creasing the sample rate. Upon decreasing the sam- 
ple rate, out of band energy is removed prior to de- 
creasing the sample rate. In this manner, out of band 
energy is prevented from aliasing into the passband 
or imaging into the stop band. 

The digital high-pass filter 240 removes the DC 
and AC power supply noise, such as 50 or 60 hertz 
components, from the incoming signal, which the 
echo canceler 140 cannot remove. 

The echo canceler 140 cancels the type of signal 
energy sent by transmitter 130 for canceling near-end 
echo. The high-pass filter removes the DC and AC 
power supply noise, such as 50 or 60 hertz compo- 
nents, prior to the received signal containing them is 
passed to the echo canceler. 

The echo canceler adapts to the impulse re- 
sponse of the signal path, including high-pass filter 
240, to eliminate the echo from the transmitted sig- 
nal. After the echo canceler has settled, its impulse 
response would look like the convolution of all of the 
combined filters in the circuit. 

Receiver 250 receives the output of the summing 
junction 260 and demodulates the signal. Receiver 
250 includes a filter 252 to remove out of band ener- 
gy. The output of receiver 250 is provided to host in- 
terface 100. 

The effect of digital high- pass filter 240 domin- 
ates the performance of modem 110. If there were 
only one filter that was a minimum phase FIR filter, 
preferably, it would be the digital high-pass filter 240. 
When all of the transmitter filter 132, digital low pass 
filter 160, digital low-pass filter 230, digital high-pass 
filter 240, and receiver filter 252 being minimum 
phase FIR filters, a two-decibel performance im- 



5 



9 



EP0 622 899 A1 



10 



provement was achieved. 

Minimum phase FIR filters provide the advan- 
tage of reduced order or length for a given magnitude 
specification. A minimum phase FIR filter has mini- 
mum length, which means both minimum time for an 5 
input signal to pass through a filter in the process of 
being filtered and minimum number of tap weight 
coefficients. A minimum phase FIR filter also has the 
advantage of lower coefficient sensitivity compared 
with the more common linear phase filters. Optimum w 
minimum phase FIR filters provide the enhanced ad- 
vantage of the shortest possible length for a given set 
of constraints. The theoretical coefficients are irra- 
tional floating point numbers. Any implementation of 
the coefficients, if finite in bit length, requires that the 15 
coefficients be represented as rational numbers. Im- 
plementation is hardware constrained to fixed-point 
arithmetic of typically 16 or 32 bits. The coefficients 
can be either positive or negative. Due to lower coef- 
ficient sensitivity, filter coefficient quantization can 20 
be achieved with less precision for an optimum mini- 
mum phase FIR filter than for a minimum phase FIR 
filter, and the filter can still be effective. 

Linear phase FIR filters have passband zeros 
some of which are on, inside, or outside the unit circle 25 
in the complex z-plane. A finite impulse response fil- 
ter is unconditionally stable regardless of where the 
zeros are located because a finite impulse response 
filter has only zeros and no poles. A strict definition 
of "minimum phase FIR filters" means there are no 30 
zeros outside the unit circle in the complex z-plane. 
A FIR filter that has all of its zeros on the unit circle 
at z=1 is a special-case filter known as a differentia- 
tor. Minimum phase FIR filters typically will have 
zeros both on and inside the unit circle. A minimum 35 
phase FIR filter is a filter in which the passband zeros 
lie on or inside the unit circle and the stop-band zeros 
lie on the unit circle. 

The following procedure is preferred for each op- 
timum minimum phase FIR filter designed. The de- 40 
sired filter frequency response is the specification to 
be met. A stop band frequency location and attenua- 
tion are specified, as are a transition band and pass 
band attenuation. These specifications provide the 
input for a commercially available software package 45 
capable of calculating the polynomial coefficients, 
along with the constraint that no zeros should be out- 
side the unit circle in the complex z-plane. The soft- 
ware should use a more general approach to deter- 
mining the linear-phase filter than the method of Her- 50 
mann and Sch ussier. Using the software package, 
the coefficients of an impulse response polynomial 
are generated. The polynomial represents the im- 
pulse response of a linear-phase filter meeting the in- 
put specifications. 55 

A soft ware package that may be used to generate 
the coefficients of an impulse response polynomial is 
Meteor. Meteor is a public-domain software package 



developed with U. S. Government funding and avail- 
able through Princeton University, Princeton, New 
Jersey. See also IEEE Transactions on Signal Proc- 
essing, August 1992, Vol. 40, No. 8, pp. 1901-1909. 

The impulse response polynomial is factored to 
obtain the roots. The roots are obtained by providing 
the coefficients of the impulse response polynomial 
as the input to another commercially available soft- 
ware package with the ability to convolve and decon- 
volve polynomials with user specified high precision. 
The polynomial factorization requires very high preci- 
sion, typically greater than 64-bit double-precision ar- 
ithmetic, to be numerically accurate. The software 
achieves the high user-specified precision by extend- 
ing precision in software rather than relying on the 
precision of the floating-point processor available to 
execute the software. 

A software package that may be used to convolve 
and deconvolve polynomials is Mathematica. Mathe- 
matics is commercially available through Wolfram 
Research Corporation, Champagne, Illinois. Mathe- 
matica enables the user to specify arbitrarily high 
precision for the root calculations. 

When a minimum phase FIR filter is being real- 
ized, the impulse response polynomial contains twice 
as many roots as are needed. The roots are zeros of 
the transfer function and occur in four possible cate- 
gories: 

(1) zeros which lie on the real axis occur in reci- 
procal pairs; 

(2) complex zeros not on the unit circle in the 
complex z-plane occur as complex conjugate 
pairs and their reciprocals; 

(3) real stop band zeros occur in doubles located 
either at z=1 or at z=-1; and 

(4) complex stopband zeros occur in double con- 
jugate pairs on the unit circle. 

The polynomial roots are either real or are com- 
plex conjugate pairs. To meet the definition of a mini- 
mum phase FIR filter, only zeros on or inside the unit 
circle of the complex z-plane are retained. The zeros 
that determine a minimum phase FIR filter are a sub- 
set of the zeros of the impulse response polynomial. 
Half of the roots are in the correct locations to form 
the minimum phase filter. The other half of the roots 
form a maximum phase filter. 

One of the double sets of zeros on the unit circle 
is retained. These zeros are the stopband zeros of the 
minimum phase FIR filter. 

The radial distance from the origin of the complex 
z-plane to a root is r k . The roots within the unit circle, 
r k < 1, are retained. These roots correspond to the 
passband zeros of the minimum phase FIRfilter. One 
root of each complex conjugate pair is located at a ra- 
dius r k from the origin of the complex z-plane, where- 
as the other root of the complex conjugate pair is lo- 
cated at 1/r k . The roots within the unit circle, r k < 1 , are 
retained for minimum phase filter realization. The re- 
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dundant roots in the passband at reciprocal locations 
outside the unit circle, r k > 1, are discarded. 

If two roots have a magnitude precisely equal to 
one, one root is retained and one root is discarded. 

After the roots to be retained are selected, a new 
polynomial, H min (z), expressed as a function of the 
complex variable z, is constructed. This new polyno- 
mial, constructed using the retained roots, is the im- 
pulse response of an optimum minimum phase FIR 
filter. The new polynomial is formed in the following 
manner by multiplying first order terms together, 
where the first order terms include the roots: 

H min (z) = (z - r^z - r 2 )(z - r 3 )...(z - rO...(z - r N ) 

H min (z) = ZN + Cl zN - 1 + C*ZN - 3 + ... + CfcZ* 
+ ... + C N . iZ 1 +c N 

the new polynomial coefficients, c 1t C2, C3...CK, are the 
tap weights or coefficients of the optimum minimum 
phase FIR filter. It is understood that complex conju- 
gate pairs may be expressed as second order terms. 
The filter thus generated is an optimum minimum 
phase FIR filter having minimum delay and minimum 
length in accordance with the minimum phase alter- 
nation theorem. 

While the degree of precision required to factor 
the original impulse response polynomial is relatively 
high, the degree of precision required to calculate the 
new polynomial is even greater. A root is a numerical 
quantity which may be expressed in any number sys- 
tem. When expressed in the binary number system, 
a root is expressed in bits. The degree of precision re- 
quired to calculate the new polynomial is a function of 
the number of roots retained and the number of bits 
in the binary representation of each root. When gen- 
erating the new polynomial by multiplying the retained 
roots together, it is desired to retain all of the informa- 
tion contained in the retained roots. For N roots, each 
having a binary representation B bits in length, the 
precision in bits, P, to retain all of the information in 
the roots is given by: 

P = 2B + log 2 N 

The product of two multiplicands having binary 
representations of bit length B results in a product 
having a binary representation (2B) bits in length. The 
sum of two such products each of length 2B bits will 
result in a number (2B)+1 bits in length. It is clear that 
each doubling of the number of roots, each having a 
binary representation B bits in length, results in one 
additional bit in the final product. 

In some applications, it may not be necessary to 
meet the minimum phase alternation theorem to ach- 
ieve an acceptable set of tap weights, with the filter 
near optimum. A nonoptimum filter can deviate from 
an optimum filter or from the minimum phase alter- 
nation theorem and still provide acceptable filter de- 
sign. Such a filter design will have more taps than a 
minimum phase FIR filter design. Factors such as the 
degree of precision of calculations, more or fewer rip- 
ples in the passband or stopband of the ideal frequen- 



cy response, and tolerance specifications of the ideal 
frequency response result in filter designs having fil- 
ter coefficients designed in accordance with the dis- 
closed method that provide some deviation from mini- 

5 mum phase FIR filters. It has been found experimen- 
tally that small deviations from the design method 
easily results in filters with two or three times as 
many coefficients. This indicates that the coefficient 
calculation method is sensitive. 

10 Figure 2 discloses a block diagram of an alternate 

embodiment modem 110' in which the converters 
may be oversampling converters. The number of sta- 
ges of interpolation and decimation will differ from 
application to application. The reference numerals 

15 with a prime designation identify elements in the al- 
ternate embodiment where the function is similar to 
the element in the embodiment of Fig. 1 without a 
prime. 

Analog-to-digital converter 220' and DAC 

20 170' operate at a rate higher than the Nyquist rate, 
such as a multiple of the Nyquist rate. For example, 
analog-to-digital converter 220* may operate at M 
times the Nyquist rate and analog-to-digital converter 
170' may operate at L times the Nyquist rate. Often 

25 M and L are equal, but the multiples are not constrain- 
ed to be equal. The oversampling analog-to-digital 
converter 220* and DAC 170' can be, but are not lim- 
ited to being noise shaping coders, sigma-delta mod- 
ulators, or multibit converters of other architectures. 

30 When oversampling analog-to-digital converter 
220' is a noise shaping coder, filter 221 is provided to 
shape the noise. Preferably, filter 221 is an optimum 
minimum phase FIR filter. When oversampling digi- 
tal-to-analog converter 170' is a noise shaping coder, 

35 filter 171 is provided to shape the noise. Preferably, 
filter 171 is an optimum minimum phase FIRfi Iter. Fil- 
ters 221 and 171 are not required when the conver- 
ters are not oversampling. 

The process of lowering the sampling rate of a 

40 signal is called decimation. Decimation filter 306 re- 
moves the noise beyond the Nyquist frequency, 
which is known as antialias filtering. Following the fil- 
tering, decimation filter 306 lowers the sampling rate 
to the Nyquist rate. Decimation filter 306 includes a 

45 wide band digital low pass antialiasing filter 361, fol- 
lowed by decimators 362 and 364 which provide one 
or more stages of decimation. 

Narrow band digital antialiasing filter 363 pro- 
vides filtering between stages of decimation. This 

so can be achieved most efficiently in multiple stages as 
disclosed in Multirate Digital Signal Processing by 
Crochiere and Rabiner, Prentice-Hall, 1983, the dis- 
closure of which is hereby incorporated by reference. 
The Crochiere and Rabiner guidelines as to the num- 

55 ber of stag es were fol lowed . 

The first stage of the decimation filter need not 
be a minimum phase FIR filter if the sampling rate in 
this stage is sufficiently high and the order of the filter 
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is sufficiently low such that the delay through the fil- 
ter is insignificant compared to the delay in the narrow 
band filters down stream in the signal path. For addi- 
tional information on sigma-delta modulators fol- 
lowed by decimators, see J. C. Candy and G. C. 5 
Temes, eds., tutorial introduction from OversampHng 
Delta-Sigma Data Converters, IEEE Press, 1992. An- 
other purpose of the decimation filter is to remove out 
of band quantization noise from the oversampling A/D 
converter 220' when A/D converter 220* is an over- 10 
sampling sigma-delta modulator or other type of noise 
shaping coder. 

Filter 210' serves a similar purpose as filter 210, 
to remove noise above the Nyquist rate. Since the 
sampling rate in A/D converter 220' is M times the Ny- 15 
quist frequency, the filter order of filter 210' can be 
lower than the filter order of filter 210, and the stop- 
band frequency of filter 210' can be higher such that 
the stopband frequency of filter 210' is M times the 
sampling frequency, less the Nyquist frequency. 20 

Interpolation filter 304 increases the sampling 
rate above the Nyquist rate and removes spectral im- 
ages that result from the interpolation process. This 
can be achieved most efficiently in multiple stages. 
Interpolation filter 304 includes narrow band digital 25 
low pass reconstruction filter 341, followed by inter- 
polators 342 and 344 which provide one or more sta- 
ges of interpolation. Wide band digital low pass recon- 
struction filter 343 provides filtering between stages 
of decimation. The last stage of interpolation filtering 30 
need not necessarily be a minimum phase FIR filter, 
given that the sampling rate is sufficiently high and 
the order of the filter is sufficiently low such that the 
delay through the last filter stage is insignificant com- 
pared with the delay through the preceding stages of 35 
narrow band interpolation filters. 

The wideband low pass reconstruction filter 
180' removes spectral images that result from the 
oversampling digital-to-analog converter 170'. If the 
oversampling digital-to-analog converter is a sigma- 40 
delta modulator or other type of noise shaping coder, 
wideband low pass reconstruction filter 180' re- 
moves out of band quantization noise. 

The filters shown in Figure 2 are distributed, hav- 
ing multiple interpolators in the transmit path 300' and 45 
multiple decimators in the receive path 302'. When 
the number of stages of decimation is I in number, the 
total distributed decimation can be represented as 

I so 

M=nMi 

i=l 

and the number of stages of interpolation is J in num- 
ber, the total distributed decimation can be represent- 55 
ed as 



L=n Lj. 

i=l 



Minimum phase FIR filters or optimum minimum 
phase FIR filters can be implemented on an integrat- 
ed circuit, such as a general- purpose or special-pur- 
pose digital signal processor. Such an implementa- 
tion would include multipliers for multiplying incoming 
discrete signals times tap weight coefficients of a fil- 
ter to produce a plurality of products. The implemen- 
tation would also include an adder or adders for sum- 
ming the products, and a register in which to tempor- 
arily store the products. 



Claims 

1. Adata converter for converting a signal from one 
form to another, comprising: 

a converter (110, 300, 302, 100', 300', 
302') for converting an incoming signal from a 
signal having a first form as it is received to a sig- 
nal having a second form; and 

at least one minimum phase finite impulse 
response filter (132, 160, 230, 240 or 252) for fil- 
tering the signal in one of said first and second 
forms. 

2. A data converter as recited in claim 1 , for convert- 
ing an analog signal into a plurality of digital sig- 
nal samples, wherein the data converter is an 
analog-to-digital converter (110, 302) for convert- 
ing an incoming analog signal into a plurality of 
digital signal samples, and the at least one mini- 
mum phase finite impulse response filter (230, 
240 or 250) is for receiving and filtering the plur- 
ality of digital signal samples produced by the 
analog-to-digital converter. 

3. Adata converter as recited in claim 2, wherein the 
minimum phase finite impulse response filter is 
a digital low pass filter (230) for lowering the sam- 
pling rate of the digital signal samples. 

4. Adata con verier as recited in claim 2, wherein the 
minimum phase finite impulse response filter is 
a high-pass filter (240) for filtering DC and power 
supply noise from the digital signal samples. 

5. Adata converter as recited in claim 2, wherein the 
minimum phase finite impulse response filter is 
a receiver filter (250) for filtering out-of band en- 
ergy from the digital signal samples. 

6. Adata converter as recited in claim 2, wherein the 
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oversampling digital-to-analog converter is a 
noise shaping coder, the noise shaping coder fur- 
ther comprising a noise shaping f i Iter (221 or 1 71 ) 
that is a minimum phase finite impulse response 
filter. 



a decimator for reducing the number of 
said first digital signal samples provided from 
said digital antialiasing filter to provide a digital 
output signal having a predetermined signal rate 
less than the rate of said first digital signal sam- 
ples; 

a digital low pass reconstruction filter for 
receiving a plurality of second digital signal sam- 
ples and for filtering images of a base band from 
said second digital signal samples; 

an interpolator for increasing the number 
of said second digital signal samples provided 
from said digital low pass reconstruction filter to 
provide a digital output signal having a predeter- 
mined signal rate greater than the rate of said 
second digital signal samples; and 

an oversampling digital-to-analog conver- 
ter (170')for converting the digital output signal 
received from the interpolator into an analog sig- 
nal for sending, wherein at least one of said filters 
is a minimum phase finite impulse response fil- 
ter. 



11. A data converter as recited in claim 7, 

or a coder-decoder as recited in claim 
10, wherein the oversampling digital-to-analog 
converter (170') is a noise shaping coder, the 
5 noise shaping coder further comprising a mini- 

mum phase finite impulse response noise shap- 
ing filter. 

12. A coder-decoder for sending and receiving sig- 
nals comprising: 

an oversampling analog-to-digital conver- 
ter (220') for converting an incoming analog sig- 
nal into a plurality of first digital signal samples; 

a digital antialiasing filter (361) for receiv- 
ing the plurality of first digital signal samples; 

a decimator (362 or 364) for reducing the 
number of said first digital signal samples provid- 
ed from said digital antialiasing filter (361 ) to pro- 
vide a digital output signal having a predeter- 
mined signal rate less than the rate of said first 
digital signal samples; 

a digital low pass reconstruction filter 
(341) for receiving a plurality of second digital sig- 
nal samples and for filtering images of the base 
band from said second digital signal samples; 

an interpolator (342 or 344) for increasing 
the number of said second digital signal samples 
provided from said digital low pass reconstruc- 
tion filter (341) to provide a digital output signal 
having a predetermined signal rate greater than 
the rate of said second digital signal samples; and 
a digital-to-analog converter (170') for 
converting the digital output signal received from 
the interpolator (342 or 344) into an analog signal 
for sending, wherein at least one of said filters is 
a minimum phase finite impulse response filter. 

13. A coder-decoder as recited in claim 10 or 12, 
wherein said coder-decoder is employed as a 

40 front-end circuit for a modem (110,110'). 

14. A coder-decoder for sending and receiving sig- 
nals comprising: 

an oversampling analog-to-digital conver- 
ts ter (220') for converting an incoming analog sig- 
nal into a plurality of first digital signal samples 
and for shifting quantization noise out of the 
baseband of said incoming signal; 

a digital antialiasing filter (361) for receiv- 
so ing the plurality of first digital signal samples and 
for filtering said quantization noise from said first 
digital signal samples; 

a decimator (362 or 364) for reducing the 
number of said first digital signal samples provid- 
55 ed from said digital antialiasing filter (361) to pro- 

vide a digital output signal having a predeter- 
mined signal rate less than the rate of said first 
digital signal samples; 



7. A data converter as recited in claim 1 , for convert- 
ing a plurality of digital signal samples to an ana- 
log signal, wherein the data converter is a digital- 
to-analog converter (170) for converting a plural- 10 
ity of digital signal samples into an analog signal 
and the at least one minimum phase finite im- 
pulse response filter (132 or 160) is for filtering 
the plurality of digital signal samples provided to 
the digital-to-analog converter. 15 

8. A data converter as recited in claim 7,wherein the 
at least one minimum phase finite impulse re- 
sponse filter is a digital low pass filter (160) that 
increases the sample rate of the digital signal 20 
samples. 

9. Adata converter as recited in claim 7,wherein the 
at least one minimum phase finite impulse re- 
sponse filter is a transmitter filter (132) for filter- 25 
ing out-of-band energy from the digital signal 
samples. 

10. A coder-decoder for sending and receiving sig- 
nals comprising: 30 

an analog-to-digital converter (220 or 
220') for converting an incoming analog signal 
into a plurality of first digital signal samples; 

a digital antialiasing filter for receiving the 
plurality of first digital signal samples; 35 
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a digital low pass reconstruction filter 
(341 ) for receiving a plurality of second digital sig- 
nal samples and for filtering images of the base 
band from said second digital signal samples; 

an interpolator (342 or 344) for increasing 5 
the number of said second digital signal samples 
provided from said digital low pass reconstruc- 
tion filter (341) to provide a digital output signal 
having a predetermined signal rate greater than 
the rate of said second digital signal samples; 

an oversampling digital-to-analog conver- 
ter (170') for converting the digital output signal 
received from the interpolator (342 or 344) into an 
analog signal for sending; and 

a digital high pass filter (240) for receiving 
the plurality of first digital signal samples from the 
digital antialiasing filter (361) and for filtering low 
frequency noise from the first digital signal sam- 
ples, wherein at least one of said filters is a mini- 
mum phase finite impulse response filter. 

15.' A modem for sending and receiving signals com- 
prising: 

an oversampling analog-to-digital conver- 
ter (220') for converting an incoming analog sig- 
nal into a plurality of first digital signal samples 
and for shifting quantization noise out of the 
baseband of said incoming signal; 

a digital antialiasing filter (361) for receiv- 
ing the plurality of first digital signal samples and 
for filtering said quantization noise from said first 
digital signal samples; 

a decimator (362 or 364) for reducing the 
number of said first digital signal samples provid- 
ed from said digital antialiasing filter to provide a 
digital output signal having a predetermined sig- 
nal rate less than the rate of said first digital signal 
samples; 

a digital low pass reconstruction filter 
(341 ) for receiving a plurality of second digital sig- 
nal samples and for filtering images of the base 
band from said second digital signal samples; 

an interpolator (342 or 344) for increasing 
the number of said second digital signal samples 
provided from said digital low pass reconstruc- 
tion filter (341) to provide a digital output signal 
having a predetermined signal rate greater than 
the rate of said second digital signal samples; 

an oversampling digitai-to- analog conver- 
ter (170') for converting the digital output signal 
received from the interpolator into an analog sig- 
nal for sending; 

a digital high pass filter (240) for receiving 
the plurality of first digital signal samples from the 
digital antialiasing filter (361) and for filtering low 55 
frequency noise from the first digital signal sam- 
ples; and 

a transmitter (130) for modulating re- 



ceived data into said plurality of second digital 
samples the transmitter, including a transmitter 
filter (132) to filter the plurality of second digital 
samples, wherein at least one of said filters is a 
minimum phase finite impulse response filter. 

16. A modem for sending and receiving signals com- 
prising: 

an oversampling analog-to-digital conver- 
ter (220') for converting an incoming analog sig- 
nal into a plurality of first digital signal samples 
and for shifting quantization noise out of the 
baseband of said incoming signal; 

a digital antialiasing filter (361) for receiv- 
ing the plurality of first digital signal samples and 
for filtering said quantization noise from said first 
digital signal samples; 

a decimator (362 or 364) for reducing the 
number of said first digital signal samples provid- 
ed from said digital antialiasing filter (361) to pro- 
vide a digital output signal having a predeter- 
mined signal rate less than the rate of said first 
digital signal samples; 

a digital low pass reconstruction filter 
(341 ) for receiving a plurality of second digital sig- 
nal samples and for filtering images of the base 
band from said second digital signal samples; 

an interpolator (342 or 344) for increasing 
the number of said second digital signal samples 
provided from said digital low pass reconstruc- 
tion filter (341) to provide a digital output signal 
having a predetermined signal rate greater than 
the rate of said second digital signal samples; 

an oversampling digital-to-analog conver- 
ter (170') for converting the digital output signal 
received from the interpolator into an analog sig- 
nal for sending; 

a digital high pass filter (240) for receiving 
the plurality of first digital signal samples from the 
digital antialiasing filter (361) and for filtering low 
frequency noise from the first digital signal sam- 
ples; and 

a receiver (250) for receiving the first dig- 
ital signal samples filtered by the digital high- 
pass filter, the receiver including a receiver filter 
(252) to filter the received first digital signal sam- 
ples, wherein at least one of said filters is a mini- 
mum phase finite impulse response filter. 

1 7. A modem for sending and receiving signals com- 
prising: 

an oversampling analog-to-digital conver- 
ter (220') for converting an incoming analog sig- 
nal into a plurality of first digital signal samples 
and for shifting quantization noise out of the 
baseband of said incoming signal; 

a digital antialiasing filter (361) for receiv- 
ing the plurality of first digital signal samples and 
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for filtering said quantization noise from said first 
digital signal samples; 

a decimator (362 or 364) for reducing the 
number of said first digital signal samples provid- 
ed from said digital antialiasing filter (361) to pro- s 
vide a digital output signal having a predeter- 
mined signal rate less than the rate of said first 
digital signal samples; 

a digital low pass reconstruction filter 
(341) for receiving a plurality of second digital sig- 10 
nal samples and for filtering images of the base 
band from said second digital signal samples; 

an interpolator (342 or 344) for increasing 
the number of said second digital signal samples 
provided from said digital low pass reconstruc- 15 
tion filter (341) to provide a digital output signal 
having a predetermined signal rate greater than 
the rate of said second digital signal samples; 

an oversampling digital-to-analog conver- 
ter (170') for converting the digital output signal 20 
received from the interpolator (342 or 344) into an 
analog signal for sending; and 

a digital high pass filter (240) for receiving 
the plurality of first digital signal samples from the 
digital antialiasing filter (361) and for filtering low 25 
frequency noise from the first digital signal sam- 
ples; 

a transmitter (130) for modulating re- 
ceived data into said plurality of second digital 
samples the transmitter, including a transmitter 30 
filter (132) to filter the plurality of second digital 
samples; 

a receiver (250) for receiving the first dig- 
ital signal samples filtered by the digital high pass 
filter, the receiver including a receiver filter (252) 35 
to filter the received first digital signal samples, 
wherein at least one of said filters is a minimum 
phase finite impulse response filter; and 

an echo canceler (140), the echo canceler 
for receiving the filtered plurality of second digital 40 
samples provided by the transmitter, and for re- 
ceiving the first digital signal samples also re- 
ceived by the receiver, said echo canceler for 
canceling from the first digital signal samples the 
type of energy sent by the transmitter whereby 45 
echo is canceled from the plurality of first digital 
signal samples, the echo canceler providing as 
an output to a summing junction the plurality of 
first digital signal samples from which the echo 
was canceled, the summing junction adding said so 
plurality of first digital signal samples to the plur- 
ality of first digital signal samples received from 
the digital high pass filter and providing the sum- 
med plurality of first digital signal samples to the 
receiver. 55 



modem as recited in claim 15,16 or 17, wherein 
the minimum phase finite impulse response filter 
is an optimum minimum phase finite impulse re- 
sponse filter. 

19. Amodem as recited in claim 17, wherein the echo 
canceler is adaptive. 



18. A data converter as recited in claim 2,6,7 or 11, 
or a coder-decoder as recited in claim 14, or a 
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